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Abstract 



We study first and second order coherence of trapped dilute Bose gases using 
appropriate correlation functions. Special attention is given to the discussion 
of second order or density correlations. Except for a small region around 
the surface of a Bose-Einstein condensate the correlations can be accurately 
described as those of a locally homogeneous gas with a spatially varying chem- 
ical potential. The degrees of first and second order coherence are therefore 
functions of temperature, chemical potential, and position. The second or- 
der correlation function is governed both by the tendency of bosonic atoms 
to cluster and by a strong repulsion at small distances due to atomic inter- 
actions. In present experiments both effects are of comparable magnitude. 
Below the critical temperature the range of the bosonic correlation is affected 
by the presence of collective quasi-particle excitations. The results of some 
recent experiments on second and third order coherence are discussed. It is 
shown that the relation between the measured quantities and the correlation 
functions is much weaker than previously assumed. 

03.75.Fi,05.30.Jp 
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I. INTRODUCTION 



The defining property of a Bose-Einstein condensate is tlie macroscopic occupation of a 
single one-particle quantum state , a condition often referred to a off-diagonal long range 
order. Even tliougfi Bose-Einstein condensation lias been long known to exist in systems 
such as superfluid ^He, the off-diagonal long range order itself was inaccessible to direct 
experimental verification. The achievement of Bose-Einstein condensation in dilute alkali 
gases [0-0], however, has revived interest in such direct experimental access. Indeed, after 
a number of theoretical studies PHTl]], it has become possible to demonstrate off-diagonal 
long range order as made manifest by an interference experiment ||T^ and that finding has 
been confirmed by excellent agreement between experiment and theory ||13|,|14|. 

The ability of a Bose-Einstein condensate to produce an interference pattern of high 
visibility is closely related to the first order coherence of an optical laser beam. More 
recently, two additional experiments have also explored certain aspects of the second ||T^ 
and third order coherence [0] of a trapped Bose condensed gas. It is therefore the aim of 
this paper to extend the theory of optical coherence fl^ to the description of trapped dilute 
Bose gases. Our main focus will be on first and second order coherence. 

The preferred tools to study coherence are correlation functions. While first order co- 
herence relies on the properties of the reduced single-particle density matrix, which are 
expressed in the first order correlation function, second order coherence depends mainly on 
density correlations. Here, the strong short range repulsion between interacting atoms adds 
interesting new aspects to the discussion of atomic coherence, not present in its optical pre- 
decessor. Density correlations are an important object of study, of course, even apart from 
their bearing on coherence. 

The paper is organized as follows. We introduce in Sec. |II A| appropriate first and sec- 
ond order correlation functions, which are then used in Sec. [II B| to define coherence. It is 
the aim of the remainder of the paper to calculate these correlation functions for an exter- 
nally trapped Bose gas. In most of our explicit calculations we assume a three-dimensional 
harmonic trapping potential and evaluate our results in the thermodynamic limit, i.e., we 
will not account for effects due to a finite atom number. This restriction is justified by the 
large number of atoms (A^ ~ 10^) in recent experiments. We thereby show that coherence 
is not only a function of temperature and chemical potential but also of the local position 
of the atoms. We consider first in Sec. |T| the case of an ideal Bose gas. This allows us 
to study some basic aspects of the problem and in particular the applicability of a local 
density approximation. Atomic interactions are added in Sec. The main modifications 
are brought about by the existence of quasi-particle like excitations and by the strong short 
range repulsion between the atoms. In Sec. [V| we discuss the coherence properties that follow 
from the correlation functions calculated earlier. We also comment on the aforementioned 
measurements of higher order coherence [0,11^ and a previous theoretical interpretation 



131. We summarize our conclusions in Sec. VI. 



II. DEFINITIONS 
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A. Correlation functions 



We describe the atoms by a quantum field ^(r, t) which obeys the equal-time commuta- 
tion relations 

[^(r,t),^^(r',t)] = 5(r-r') (2.1) 
[^{r,t),^{r',t)] = 0. (2.2) 

In general, physical properties of the gas can be expressed in terms of correlation func- 
tions, i.e., expectation values of the field operators of the type 

{^\r,,h)^\r2, h)--- ^(r'2, 4)^(r;, t\)). (2.3) 

Here, a description of the dynamics using the Heisenberg picture has been tacitly assumed. 
However, we shall restrict ourselves to discussing spatial correlations at thermal equilibrium 
and at a fixed time. We therefore omit further mention of the time arguments. Temporal 
correlations have been discussed recently in the context of condensate phase fiuctuations 

]ig,[n]jg-g5i. 

For later convenience we introduce a special notation for the correlation functions 

G(^)(r,r') = (^^(r)^(r')) (2.4) 
G(2)(r^r') = (^t(r)^t(r')^(r')^(r)). (2.5) 

The first order correlation function G^^^ (r, r') contains the same information as the reduced 
one-particle density matrix which furnishes the expectation values of one-particle observ- 
ables like the momentum or position density. The formal relation between first order corre- 
lation function and reduced one-particle density matrix is given by G'^^\y, r') = p*^^^(r', r) = 
(j'\p^^^\y). The first order correlation function can alternatively be expressed in terms of the 
Wigner function 

W^'P^^) = j d^^-'"''"' (2.6) 

which is a quantum mechanical analog of a classical single-particle phase space distribu- 
tion. The connection between the quantum mechanical correlation function and statistical 
mechanics is indeed illustrated by the Fourier transform relation 

G«(r,r')= /rfpe-'P('-'-')/^l^(p,^^^). (2.7) 



In contrast, the second order correlation function G'-^''(r, r') expresses the joint probability 
of detecting two different atoms at the locations r and r' respectively. Its connection with 
the density correlation function is given by 

{i>\v)i!{Y)¥{v')i>{Y')) = G^^\t, y') + G(^)(r, r) 5{y - r'). (2.8) 

It is because the function G^^^ is defined as the average of a normally ordered product of field 
operators, that it expresses the correlation of distinct pairs of atoms. The delta function 
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contribution in Eq. ( p.8| ) represents by contrast the autocorrelation of individual atoms. 
Discussions of the coherence properties of fields are most conveniently phrased in terms 
of factorization of their correlation functions WTl. The autocorrelation contribution in Eq. 



however stands in the way of any rigorous factorization of the density correlation into 
a product of functions of r and r'. We shall therefore concentrate on discussing the normally 
ordered correlation function G'*^^^(r,r') which can indeed meet the factorization condition 
under suitable physical circumstances. 



B. Coherence and Bose-Einstein condensation 

In the optical context coherence refers to a state of the system with interference properties 
that are as close as possible to those of a precisely defined classical field. In more technical 
terms, successive orders of coherence can be described by means of factorization conditions 
imposed upon a succession of field correlation functions. The analog of n-th order coherence 
for the atomic field would require all normally ordered correlation functions of order m < n 
to factorize into a form 

{4^\r,)4^\T2) ■ ■ ■ ^{r',)^{r[)) = V*(ri)z^*(r2) ■ ■ ■ ^{r',)^{r[), (2.9) 

i.e., into products of the same complex valued function ipi'"^)- If fi^ld were in a pure 
coherent state, the factors ipi'^) could be regarded as field-expectation values (\l/(r)). This 
identification, however, is by no means necessary to secure the factorizations of Eq. ( p. 91) . 
It would even contradict the strict conservation of the total atom number. 

To identify the function ip{r) we begin by considering first order coherence, which secures 
the maximum effects of intensity interference between spatially separated parts of the system. 
In this sense, perfect first order coherence 

^^(r,^ =^*(r)^(r') (2.10) 

with 

^(r) = v^(r|</.) (2.11) 

will hold if all atoms occupy the same one-particle state 10). Here denotes the total atom 
number of the system. This definition is an idealization of the Penrose and Onsager definition 
of Bose-Einstein condensation which requires the reduced one-particle density matrix to 
have a single macroscopic eigenvalue. The function ip{r) can thus be introduced without 
breaking the gauge symmetry associated with the conservation of total atom number. The 
modulus of ipi^) can be reexpressed in terms of the first order correlation function as 

|V.(r)| = v'GW(r,r). (2.12) 

In addition to the factorization of G^^\ perfect second order coherence implies the factor- 
ization condition 

G(2)(r,rO = |^(r)n^(r')P (2.13) 
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so that 



G^^\r,r') = G'^^\r,r)G^^\r',r'). (2.14) 

First order coherence does not, of course, require second order coherence to hold. It does, 
however, imply a weaker factorization property of the function G^'^^ which includes the 
statement 

G^'\T,v')=g2mr)\'mv')\^ (2.15) 

where g2 is a non-negative real number. For example, a first order coherent system with 
definite total atom number can achieve second order coherence, i.e. g2 = 1, only to an 
accuracy of order 0{1/N). 

In fact, any real atomic system can exhibit coherence only in an approximate way. In 
order to define a local measure of coherence we introduce correlation functions which are 
normalized to attain unit modulus in the case of perfect coherence. The degree of first order 
coherence, in this sense, is expressed by 

^'■'(^■o^ (2.16) 

while 



GW(r,r)JG(i)(r',rO 



^ ^'^ G{i)(r,r)G'(i)(r',rO ^ 

is a certain measure of second order coherence. In the language of a classical field theory 
g^^> characterizes local fluctuations of the phase of the complex field amplitude, while 
related to fluctuations of its modulus. The most important property of ^7(1) is its relation to 
the contrast achievable in an interference experiment [^. The function g^'^\ on the other 



hand, expresses the tendency of atoms either to cluster or to remain spatially separated. 

In many experimental contexts, it will be difficult, to collect enough data to determine 
the dependence of the correlation functions on the two coordinates r and r'. It may then be 
simpler to measure volume integrated correlation functions 

G«(r) = |dRG«(R-^,R+^) (2.18) 

G(2)(r) = |rfRG(^)(R-^,R+^), (2.19) 

in which a spatial integration over the mean coordinate has been performed. In condensed 
matter physics, for example, scattering experiments are often used to obtain information 
about the volume integrated second order correlation function. We define normalized ver- 
sions of these functions 

s»>(r) = /.RGa.(R-.R+|) ^^^^^ 



JdR ^/G(i)(R - I, R - I) ,/G(i)(R + |, R + f ) 



r ^ 



g(2)M = /rfR G'(^)(R-|,R^^; 



/rfRG«(R - f , R - f ) ^^(R + |, R + I) ' 
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which can also be interpreted as weighted volume averages of the normalized correlation 
functions of Eqs. (|2.16| ) and ( p.l7|) . 

It will be the aim of the following sections to calculate the correlation functions we have 
defined for trapped dilute Bose gases at finite temperatures. The main emphasis will be put 
on systems that are large enough to allow a local density approximation. 



III. IDEAL GASES 

The basic concepts for calculating the correlation functions of a trapped gas are most 
easily described for ideal gases with temperatures above the critical temperature. The 
presence of a Bose-Einstein condensate at temperatures below the critical temperature gives 
rise to additional features which will be discussed in Sec. 1111 B|. The influence of atomic 



interactions on both uncondensed and condensed gases will be later dealt with in Sec. IV 



A. Temperatures above Tc 

In the simplest case, we consider an ideal Bose gas at temperature T trapped by an exter- 
nal potential ^(r). In the grand canonical ensemble different one-particle energy eigenstates 
Mn(r) of the trapping potential are populated independently. As a result we can write 

G«(r,rO =E<W«m(r')(aLam) (3.1) 

m 

G(^)(r, r')= ^ <(r)Mi'(r')Mm(r') ^in(r) (ak^i^mfln) 

klmn 

= (^t(r)^(r))(^t(r')^(r')) + | (^t(r')^(r)) p (3.2) 



The statistical independence we have assumed for the mode occupations would not, on the 
other hand, hold precisely for a canonical ensemble, i.e., a system with a fixed total atom 
number. However, the error in the correlation function induced thereby is only of the order 
0{1/N) and can be neglected for large systems. 

The occupation fluctuations characteristic of the grand canonical ensemble make the last 



line of Eq. (3^) vanish. The second order correlation functions can then be written as 

G(2)(r, r') = ^^(r, r) ^^(r', r') + |G'W(r, r')P (3.3) 
^?(2)(r,r') = l + |^?«(r,r')p. (3.4) 

The second term of G*^^^ or g'^'^'> is a quantum statistical exchange term. It manifests a 
tendency of bosonic atoms to cluster together. In the case of fermions this term would have 
a negative sign reflecting the influence of the Pauli exclusion principle. 

Since we may freely use the grand canonical ensemble at temperatures above Tc, we need 
only determine the first order correlation functions G*^^-'(r, r') and ^''^^■'(r, r'). Furthermore, no 
condensate exists in this regime and the normalized correlation functions have the limiting 
values 
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lim g^^\r,r' 







r— r' ^oo 



lim g^'^\r, r') = 1. 

|r— r'|— >oo 



(3.5) 
(3.6) 



In the following we compute the length scale of r — r' over which the correlation functions 
decrease. 



1. Three-dimensional harmonic oscillator potential 

As a first step we calculate the first order correlation function for the case of a spherical 
harmonic oscillator potential ^(r) = ma;^r^/2 with angular trap frequency uj. We therefore 
occupy the harmonic oscillator states in Eq. ( p.l|) with a Bose-Einstein distribution, so that 



GW(r,rO = E«m(r)u^(r')r^ 



(3.7) 



with the inverse temperature (3 = l/Zc^T, the fugacity z = e^'^, and the chemical potential 
fi. The energies e^n have the zero point energy, 3huj/2, removed in order to limit the fugacity 
to values between zero and one, as in a homogeneous gas. If the summand of Eq. ( p.7| ) is 
expanded in powers of z, the remaining sum over the state indices m can be compactly ex- 
pressed in terms of the single-particle propagator Go{r', r, t) = (r'| exp{—iHt)\r) for motion 
governed by the harmonic oscillator Hamiltonian, H. We obtain 



G«(r,r') = E<(r; 



Ury,[r ze 



-/3en 



m 

oo 



E 

k=0 



exp{3k(3tuj /2) Go{r', r, t = -thkf3) 

k=l 

oo 

Goir',r, kphu) 



(3.8) 



k=l 



where 

Go(r',r,r) 



muj 



3/2 



exp 



muj (r^ + r'^)(coshr — 1) + (r — r')^ 
h 2 sinh r 



(3.9) 



is derived from the Green's function for the three-dimensional harmonic oscillator p8|. The 



fugacity z must be determined self-consistently by using the relation p9| 



^=^7 

ei 1- 



k=l 



(3.10) 



in order to insure the proper total atom number N. 

A considerable simplification can be achieved if the energy spacing of the trap levels is 
much smaller than the thermal energy, i.e., f3huj -C 1. In that case Gq reduces to 



G'o(r', r, r) 



muj 



3/2 



exp 



moo (r^ + r'2)r/2 + (r - r') V'^ 
~J 2 



(3.11) 
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and the sum for G^^^ can be written in the form 



1 ^ 



E 

k=l 



exp 



A: 



3 93/2 exp 



[V-(r) + ViT')]/2 
^-[V{v) + V{v')]/2 



exp 



r — r 



/\2 



— TT- 



A 



i/fc 



knT 







,exp 


Aj^ 



where we have used the thermal wavelength 



At = h \J 27r / mkBT . 
In addition we have introduced the generalized Bose function 

oo k 1/k 
I \ \ ^ y 



fc=i 



(3.12) 



(3.13) 



(3.14) 



which is related to the commonly used Bose function (?a(x) = Y^^=\^^ jkP" and to the Rie- 
mann zeta function ^(a) = Y^'k=\ 1/^" through 



9cAx) = ga{x, 1), C(a) = 9a{i) = 9a{i, !)• 



(3.15) 



2. Local density approximation 



The implications of Eq. ( |3.12| ) can be made clearer by deriving an approximate form 
for it from Eq. (|2.7|) . There G^^^ was expressed in terms of the Wigner function. We now 
employ the local density approximation which assumes the Wigner function, W{p,q), to 
be locally identical to the momentum distribution of a spatially homogeneous system with 
constant potential energy equal to the local value V{q) 



W{p,q) 



{2nhf exp [(p2/2m + F(q) - ^i)/kBT] 



(3.16) 



This approximation implies that the atoms see an effective chemical potential /i(q) = /i — 
V{(\) which varies with position. If we Fourier transform this Wigner function as required 
by Eq. (|2.7|) the correlation function can be expressed as 



1 



G(^)(r,r') = ^^3/2 (exp 



^l-v[{Y + v')/2] 

knT 



exp 



-7r- 



Aji 



(3.17) 



The global chemical potential must then be chosen to secure the correct total atom number 

N= /drGW(r,r). (3.18) 



We can now compare Eq. ( ^.17] ) with Eq. ( |3.12| ). These expressions for G^^-* in fact 
coincide for equal arguments r = r'. However, for r 7^ r' they treat the potential energy in 
a slightly different way. The accuracy of the local density approximation evidently requires 
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that the external potential vary slowly on a length scale given by Xt- For a relative separation 
r — r' of this magnitude, the difference of the potential energy between Eqs. ( ^121) and (FTTD 
is only of the order {hu)'^ /ksT and can thus be neglected. This result can be no surprise if 
we consider that the thermal wavelength, written as 

than the natural length unit Oq = 



2TThu/kBT ao, becomes smaller 

Ifi/muj of the trapping potential precisely in the limit 
huj <^ ksT. For the trap parameters of Ref. |12], the thermal wavelength assumes a value 



of At = 0.4 /im at the critical temperature, which is three to ten times smaller than the 
harmonic oscillator length units of the anisotropic trapping potential. 

We have thus shown that both the density distribution and correlations implicit in G'^^^ 
are well described by the local density approximation, as long as the separation of the trap 
levels is small compared to the thermal energy. Of course, application of the local density 
approximation is not restricted to an isotropic harmonic oscillator potential. It can be 
applied to any trapping potential V^(r) as long as the thermal wavelength is small compared 
to the length scale of the spatial variations of the potential. Since the correlation function 
decreases on the scale of the thermal wavelength it is consistent to use for small relative 
separations (r < oq) the alternative definitions 



JdR G«(R 



|,R+ I) 



JdRGW{K, R) 

/rfRG(2)(R- |,R+ 1; 

/rfR|G(i)(R, R)|2 



(3.19) 
(3.20) 



for the normalization of the volume integrated correlation functions. 

For a harmonic oscillator potential the integration over the mean coordinate R = (r + 
r')/2 of Eq. ( |3.17| ) can easily be carried out by using the series expansion of (73/2. The 
results for the volume integrated correlation functions of Eqs. ( |2.18|) and ( p.l9| ) are then 



flUJ 



93 exp 



/i 



knT 



,exp 



— TT 



Aj' 



93 (exp 


[kBT\ 


,exp 


"-vr ""^ ' 


) 


93 (exp 


\ 

[kBT 


) 



(3.21) 
(3.22) 



The latter formula applies equally well to the case of an anisotropic harmonic potential when 
the replacement to = {cUxLUyUJzy^^ is made. As in the well known case of a spatially homo- 
geneous system, the correlations decay on a length scale given by the thermal wavelength 
At. However, the Bose function appropriate to a harmonic oscillator potential is closer 
to a linear function of its argument than the function (73/2 characteristic of a homogeneous 
system. The properties of a Bose-Einstein distribution are thus less visible in the volume 
integrated correlation function of a harmonically trapped system than in that of a homoge- 
neous system, since the atoms away from the trap center experience a lower local chemical 
potential. In view of this reduced tendency towards condensation for larger radii it may 
be an acceptable approximation to assume a Maxwell-Boltzmann occupation distribution of 
the trap levels, or equivalently to restrict the summation in Eq. ( p.l2| ) to its lowest order 
term. In this approximation we obtain the simple expressions 
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^2 \ / 



g^'\r) = exp (^"^^j , o'^'K^) = 1 + exp (^-2vr^ J . (3.23) 

Numerical calculations confirm (cf. Fig. |1|) the accuracy of these simple expressions even at 
temperatures very close to the critical temperature. 



B. Temperatures below Tc 
In Sec. 1111 A| we have calculated the first and second order correlation functions for 



uncondensed systems. The presence of a Bose-Einstein condensate, however, will lead to 
significant changes, especially in the second order correlation function. 

In principle, the presence of a Bose-Einstein condensate has no effect on the previously 
discussed expression, Eq. (|3.1|), for the first order correlation function G^^\ However, due 
to the relatively long ranged coherence of the condensate we have to perform the local 
density approximation in a slightly different way. We therefore introduce a Wigner function, 
Wt{p, q), which describes only the uncondensed atoms. The condensate, which corresponds 
to the ground state of the bare external potential, is still considered to have a wave function 
ip{r) that obeys the Schrodinger equation 

/ Ji^V^ \ 

V(r)-/iU(r) = 0, (3.24) 



y 2m 

with ^ = h{ijJx + ujy + ujz)/2 corresponding to the zero point energy of the oscillator poten- 
tial. The first order correlation function can then be decomposed into a condensed and an 
uncondensed or thermal part 

G(^)(r,r') =^*(r)^(r') + G?(r,r') (3.25) 

with 



G(!)(r,r') = y'rfpe-^P-('-'-')/^WV(p,^^). (3.26) 



Again, the thermal Wigner function Wt describes the uncondensed part of the full Wigner 
function. The thermal correlation function G^'* is given by Eq. ( |3.17| ), since the momentum 
integration which led to Eq. ( p.l7|) does not account for the existence of a condensate. 



However, in contrast to Eq. ( ^.24| ), the chemical potential in Wt is set equal to zero, (/i = 0), 
as it must be in the local density approximation. The occupation of the condensate mode, 
furthermore, has to be determined from the normalization condition 

N = jdv\^{v)\^ + JdrG?{r,r). (3.27) 

The limiting values of the normalized first order correlation function 

(7«(r,r) = l (3.28) 
lim ,W(r,rO= , ^^^^^^^'^^ (3-29) 



10 



differ from those of an uncondensed gas (cf. Eq. (|3.5| )) by the possibihty of off-diagonal long 
range order. The existence of a condensate causes the correlation length of g^^^ to be much 
larger than the thermal wavelength At- In the case of a spatially homogeneous gas the 
correlation extends throughout the entire system. The degree of first order coherence for 
large separations is then just given by the local condensate fraction. 

The second order correlation function is affected by Bose-Einstein condensation in a 
more interesting way, as can be seen in Eq. ( |3.2| ). The fluctuations of the number of atoms 
in the condensate mode depend sensitively on the statistical ensemble that is assumed for 
the calculation. The fluctuations characteristic of the grand canonical ensemble greatly 
overestimate the condensate fluctuations in any experiment with a fixed total atom number. 



Instead we assume canonical condensate fluctuations which are of the order 31 



(dlalaoao) = (aoao)(ao^o) + 0{N). (3.30) 

In the limit of a large system we can neglect the last term since it is much smaller than the 
first term which is of the order 0{N'^). To calculate the second order correlation function 
we refer back to Eq. ( p^) and note that the fluctuations of the uncondensed modes can be 
well approximated by using the grand canonical ensemble. For the condensate mode, on the 



other hand, we can write Tp{r) = \J {olaQ) Uo{r), so that we have 

G'(2)(r,r') = G(^)(r,r)G«(r',r') + |G«(r,r')|' - |^(r)|2 |^(r')|' (3.31) 

The limit for large relative separations is again given by 

lim g^^\r,r') = l. (3.33) 

r— r'|— »oo 

A certain measure of second order coherence is 

This function, which is determined by the ratio of the local condensate density to the total 
density must reduce to unity everywhere in order to describe perfect second order coherence. 
In general, it can only be close to one near the center of the trap, where the condensate 
density is large. The coherence properties of the thermal cloud outside the condensate, i.e. 
where ip{r) = 0, are not greatly different from those of an uncondensed gas. 



IV. STRONGLY INTERACTING, DILUTE GASES 

In Sec. m we have studied the basic properties of the first and second order correlation 
functions of a trapped ideal gas. In this section we discuss the way in which atomic inter- 
actions modify these results. In order to deal with recent Bose condensation experiments 
we focus on gases which are very dilute, but which contain atoms with a strong mu- 
tual repulsion. In this section we discuss only Bose condensed gases. The properties of an 
uncondensed gas can easily be obtained as a specialisation of the following treatment. 
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A. First order correlation function 



In the presence of atomic interactions, the condensate wave function obeys the finite 
temperature Gross-Pitaevskii equation (cf. Eq. ( p.24| )) 

+ V{v) + U[\iP{r)f + 2nTir)]- fx]ij{r) =0 (4.1) 



y 2m 

where the atomic interaction strength is measured hj U = A-Kfi^a/m^ with a denoting the 
two-body scattering length. The local density of thermal atoms nr(r) is given by the still- 
to-be determined thermal part of the first order correlation function (cf. Eqs. (|3.25| ) and 

(B)), 

nT(r) = G'y)(r,r) = j cipWrM. (4-2) 

In the case of an interacting gas, the local density approximation has to be applied not only 
to the distribution of thermal atoms but also to the condensate wave function in order to 
be consistent. Otherwise, the local kinetic energy of the condensate wave function could 
be larger than the kinetic energy of the low lying elementary excitations. Details of this 
argument will be presented later. We have therefore to neglect the kinetic energy term in the 
Gross-Pitaevskii equation ( [4.1D , and that is appropriate if the kinetic energy is much smaller 
than the interaction energy. Under present experimental conditions this condition is 
well fullfilled in the core of the condensate. However, the local density approximation can 
not be applied in a small region around the surface of the condensate, where the interaction 
energy ceases to dominate the kinetic energy. The radial extent of this region is usually 
called the healing length. It has a fixed thickness and therefore a vanishing relevance in the 
thermodynamic limit. This local density or Thomas- Fermi approximation has already been 
used frequently for the description of Bose-Einstein condensation experiments. In this limit, 
the condensate wave function is given by 



for positive arguments of the square root and by zero otherwise . 



Since the local density approximation fails close to the surface of the condensate we are 
only able to study correlation functions with both of their spatial arguments located within 
the condensate or outside of it. Inside the condensate, the first order correlation function is 
given by 

G«(r,r') =^*(r)V^(r') + Jdpe-'^-'-''-'''^/^WT{p,'^^) (4.4) 



with the condensate wave function defined by Eq. (|4.3|) . For the calculation of the thermal 
Wigner function we have to take into account that the low lying excited states of the inter- 
acting Bose gas are excitations of collective quasi-particle modes rather than single atoms. 
To be more specific, a proper treatment of the atomic interactions in the presence of a Bose- 
Einstein condensate starts with decomposing the atomic field operator into a condensate 
and a fluctuation part 
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^(r) = ^(r) + 5^(r) 



(4.5) 



with ip{r) obeying the Gross-Pitaevskii equation ( [4.1| ). Due to this choice of ^/'(r) only second 
order terms in the fluctuations 5\i/(r) and higher contribute to the interaction Hamiltonian. 
Within the Popov approximation third and fourth order terms in 5\l/(r) are treated in an 
averaged way, neglecting anomalous expectation values like (5^(r)5\l/(r)). This extends the 
idea of the Bogoliubov approximation to finite temperatures since third and fourth order 
fluctuation terms are completely neglected in the latter. It is then possible to introduce 
quasi-particle operators aj which diagonalize the approximate Hamilton operator of the 
interacting system to 



H 



(4.6) 



The quasi-particle operators aj obey Bose commutation relations. In a spatially homoge- 
neous system with volume V they are related to the fluctuation operator (5^(r) by |]33 



(5^(r) 



1 



ip-r/h 



(4.7) 



By virtue of the local density approximation both the coefficients u{p, r) and f (p, r) as well 
as the Hartree-Fock energy 



6hf(p, '■) = ^ + ^ W + 2 U Mr)]' + nrir)] - 



(4.8) 



and the Popov energy 



e(p,r) 



{umr)\ 



21 2 



2m ~'~ 



1 2 



-[u 



(4.9) 



become functions of the position r. The local spectrum of elementary excitations e(p, r) 
coincides with the well known Bogoliubov spectrum for a Bose condensed homogeneous gas 
close to zero temperature. The external potential and the interactions between uncondensed 
atoms influence the excitation spectrum only indirectly through the local condensate density 
(cf. Eq. ( |4.3| )). If the kinetic energy of an elementary excitation with momentum p is smaller 
than the interaction energy of the condensate atoms, these excitations consist of phonon- 
like quasi-particles. For large kinetic energies, the Popov energy e(p, r) coincides with the 
Hartree-Fock energy ejjp(p,r) and the elementary excitations then consist of real particles 
instead. 

Eq. (O) also explains why it was necessary to extend the local density approximation 
to the condensate wave function. Had we not dropped the kinetic energy of the condensate 
from the definition of the chemical potential in Eq. ( [4.1| ), the Popov energy would be given 
by 
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\ 



2m 



+ 



^*(r)V2^(r) 
2^ WyW 



[t/|^(r)P 



(4.10) 



For a positive kinetic energy of the condensate wave function the argument of the square 
root could become negative, leaving the Popov energy undefined. One might try intuitively 
to remedy this problem, for example, by restricting the phase space available to the quasi- 
particles. However, numerical comparisons of several such schemes have yielded even less 
agreement with the results of an exact Path Integral Monte Carlo calculations [^| than 
the complete neglect of the quasi-particle excitations that corresponds to the Hartree-Fock 
method. 

Finally, we obtain the thermal Wigner function by assuming the quasi-particle states to 
be occupied by a Bose-Einstein distribution and by expressing the mode functions u(p, r) 



and f (p, r) in terms of energies |30| 



WT(p,r) 



1 



1 



u (p,rj 
eHF(P'r) 



V (p,rj 



exp[e(p, v)/kBT] - 1 



{2'KhY e(p,r) exp[e(p, r)/A;Br] 



(4.11) 



Outside the condensate, collective effects do not play any role in determining the spec- 
trum of elementary excitations. The first order correlation function is then given approxi- 
mately by the formula ( p.l7| ) for an uncondensed ideal gas 



r, r 



^^73/2 (^exp 



M[(i- + rO/2] 



exp 



,>'^2 



-TT- 



(4.12) 



in which the local chemical potential /i(r) = /i — V{v) — 2 UnT{r) is reduced by the mean- 
field interaction between the uncondensed atoms. However, these interactions are usually 
very weak and the correlation function is then essentially that discussed in Sec. |T|. 

A numerical plot of the volume integrated, normalized correlation function g^^\r) for 
an interacting Bose condensed gas is given in Fig. The effective strength of the atomic 
interaction is best expressed in terms of the universal scaling parameter rj = fiT=o/kBTj^ 
30| , where T° is the critical temperature of the ideal gas. This scaling parameter assumes 



values between 0.3 and 0.4 in most recent experiments. As in the case of a condensed ideal 
gas, the correlation function decreases for relative separations much larger than the thermal 
wavelength, but smaller than the total size of the condensate, to the overall condensate 
fraction, i.e., the system posesses an off-diagonal long range order. However, for an inter- 
acting gas, the tail of the correlation function decreases much more slowly, as a result of 
the presence of the long wavelength quasi-particle excitations. The unintegrated correlation 
function G(^^(r, r') depends strongly on its mean coordinate R = (r + r')/2 since the local 
condensate fraction varies significantly over the atomic sample. 



B. Second order correlation function 

Atomic interactions affect the second order correlation function in a more fundamental 
way than just by changing the spectrum of elementary excitations. Atoms repel each other 
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strongly at small separations. It is therefore impossible to find two atoms at exactly the 
same location, i.e., G*^^'*(r, r) = 0. This effect is clearly not accounted for by Eqs. (p. 31 ) 



and ( |3.32| ). In order to study the influence of this hard core repulsion we briefly revisit the 



derivation of the time-dependent version of the Gross-Pitaevskii equation (|4.1|) . 

The condensate dyad ^jJ{r)^p*{r') is part of the reduced one-particle density matrix p^^\ 
which contains the same information as the first order correlation function. The dynamics 
of p^^-* depends on the reduced two-particle density matrix p^"^^ through the Schrodinger 
equation 

|p^^^ = -^[^o,i,p(^)]-^tr,{[if„p(^)]} (4.13) 

with the one-particle Hamiltonian ifo,i and the true atomic interaction potential U{ri — r2) 
in the coordinate representation 

Ho,i = -—^ + Vir,), Hi = U{vr-V2). (4.14) 
Im 

The dynamics of the second order correlation function is similary linked to the third-order 
correlation function and so on, the entire set of N coupled equations constituting the so- 
called BBGKY hierarchy. If it is possible, however, to express the reduced two-particle 
density matrix p*-^-* in terms of the reduced one-particle density matrix p^^\ the hierarchy 
can be truncated by turning Eq. (|4.13| ) into a closed nonlinear equation. 

In most treatments, the time-dependent analog of the Gross-Pitaevskii equation ( [4.1| ) 
is obtained by replacing the true atomic interaction potential by the point interaction 
Hi = U 5{yi — Y2) and by factorizing the reduced two-particle density matrix analogously 
to Eq. (|3.31|) . Even though such an application of the first order Born approximation to 
a suitably chosen pseudopotential properly reproduces the Gross-Pitaevskii equation, its 
implications for quantities other than the reduced one-particle density matrix require more 
careful examination. For example, the above treatment implies that the second order corre- 
lation function 

Gf^ (r, r') = G(^) (r, r) G^^^ (r', r') + \G^^^ (r, r') |' - |V^(r) ^ |V^(r') |' (4. 15) 

corresponds to an ideal gas of quasi-particles with G'*^^)(r, r') given by Eq. (^.4]). Such a con- 
clusion, however, incorrectly ignores the strong repulsion between atoms at small distances. 



We therefore outline a more detailed derivation of the Gross-Pitaevskii equation ||23|j35 
which also describes the two-particle density matrix appropriately. For this purpose we 
decompose the two-particle density matrix, p*^^^ = p^^^ + ^p*-^-*, into a part pg^'' representing 
products of the reduced single-particle density matrix (cf. Eq. ( |4.15| )) and a part 5p'^'^\ which 
corrects for the true two-particle correlations. Rather than neglecting the correlation part 
5p''^^ altogether, we turn Eq. ( [4.13| ) into a closed equation by solving the dynamics of Sp^"^^ 
analytically and by reexpressing the correlations, 5p'-^\ in terms of Pq'\ Of course, certain 
approximations have to be made in order to achieve a convenient analytical expression for 
5p(2). 

The equation of motion of 5p*^^-* is simplified by the assumption of an atom density 
low enough to permit neglect of three-body and higher order collisions. The two-particle 
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correlations are therefore the result of two-particle scattering only. In the context of many- 
particle physics this approach is often called a ladder approximation since the two-body 
collision is calculated up to arbitrary orders in the true interaction Hamiltonian Hj. It is 
worth noting that a literal point-like potential does not lead to any two-particle scattering, 
and that is consistent with the assumption p*^^^ = Pq'\ but it is not at all realistic. 
If we describe the two-particle scattering by the M0ller operator 



n = hm exp[-2(^ro,i + ^^0,2 + Hi)t/h] exp[i{Ho,i + Ho,2)t/h], 



(4.16) 



it can be shown that the properly approximated two-particle density matrix is given by 



P 



(2) 



Q p^^ . Thus the equation of motion for p^^^ reads | 



^pW 



^tr2 fip(')l]t]} 



(4.17) 



The M0ller operator Q is related to the two-body T-matrix by the relation T = Hj Q. The 
matrix elements of the two-body T-matrix in momentum representation are identical to 
the scattering amplitude. It is therefore the T-matrix rather than the atomic interaction 
potential Hi which should be approximated by a delta function T = U 6{ri — r2). To lowest 
order in the scattering length, this equation of motion can be further reduced to 



d_ 
di 



P 



(1) 



h 



[Ho,i,P 



(1)1 



tr2{[T,p?)]} 



(4.18) 



which coincides with the previously outlined mean-field approximation. The neglected terms, 
which are quadratic in the scattering length, describe small thermal energy shifts and the 



dissipative effects of collisions such as thermalization [23 



We conclude from the above derivation of the Gross-Pitaevskii equation that the second 
order correlation function that properly accounts for hard core repulsion is given by 



(4.19) 



The correlation function is here defined as the operator G*-^-* = p^"^^ ^ with the understanding 
that its coordinate representation corresponds to G(^)(r,r'). 

For reasons of simplicity we calculate G*^^) first for an uncondensed gas. The scattering 
involves only the wave function for the relative coordinate r = ri — r2 of the two atoms. If 
the incoming states are plane waves the effect of the M0ller operator is given by 



g«{pi ri+p2-r2)/?i _ giP-R/7i 



Jp-r/h 



(r > a) 



(4.20) 



with the center-of-mass coordinate R = (ri + r2)/2, the total momentum P = pi + p2, 
and the relative momentum p = (pi — P2)/2. Here, we have assumed that all relative 
momenta are small enough to permit taking only the momentum independent part, —a, of 
the scattering amplitude into account. With this approximation we have neglected terms of 
order pa/h which is justified if a <^ At. In a harmonic trap the ratio a/Ay is given by 



4 1 



At 15v^C(3) 



V 



5/2 



1/2 



(4.21) 
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It is mainly determined by the scaling parameter t] = fiT=o/kBTc [0 which, as we have 
noted, assumes values between 0.3 and 0.4 in most recent experiments. As a consequence, 
we can generally assume a ratio of a/Aj- ~ 0.01 at the critical temperature. 

In agreement with the local density approximation, we evaluate G^^-* for a spatially 
homogeneous gas with volume V. Denoting the relative coordinate once more as r and the 
relative momentum as p, we obtain for an uncondensed system (cf. Eq. (|3.3|) ) 



Pl,P2 



1 



Pl,P2 



E e'P"/^-a^ 



^ipr/h \ t 



Jp-r/h 



— a- 



Pl:P2 



1 

y2 



Pl:P2 \ 



^ipr/h \ t 



oipr/h \ 

r J 

pipr/h 



r 



(^Pl)(^P2) 
(^Pl)(^P2)- 



(4.22) 



After multiplying this expression out term by term and making use of the spatial isotropy 
of the homogeneous gas we can write the second order correlation function as 



G(2)(r,r') = G«(r,r) G«(r',rO (^1 + j^^^ + \G^'\r,r')\' (^1 

/ \ 2 

= G«(r,r)G«(r,r') + |G«(r,r') 
+ [G«(r,r)GW(r',r') - |G«(r,r') 



4a 



r — r 



r — r 



- 1 



r — r 



(4.23) 



an expression which should also be valid for a sufficiently slowly varying inhomogeneous 
system. The first line of Eq. (|4.23|) multiplies the correlation function of an ideal gas by the 
squared relative wave function in the limit of zero momentum scattering. The additional 
term on the second line accounts for the fact that the atoms have finite relative momenta 
with a distribution implicitly characterized by ^'•^•'(r, r'). This correction term can only be 
neglected for very small a/Ar ratios. Furthermore, only linear terms in a/|r — r'| are relevant 
for distances |r — r'| ^ a. 

The extension of the previous calculation to Bose condensed gases is straightforward 
though it requires a little more algebra due to the existence of quasi-particle excitations. 
However, most of the required modifications are absorbed into the definition of the first 
order correlation function G^^\ In addition one has to take into account the fact that the 
second order correlation function of a condensed ideal gas is given by Eq. ( |3.31|) rather than 
by Eq. ( p.3[) . With these modifications the second order correlation function of a condensed 
interacting gas is given by 



r, r 



G, 



(2), 



r, r 



1 



(4.24) 
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G'W(r,r) G'W(r',r') - \G^'\r,r') 



r — r 



r — r 



1 - \Q^^Hr,r') 



r — r 



(4.25) 



In Fig. we compare the volume integrated normalized correlation functions g^'^'' and g^Q^ 
for an uncondensed system at the critical temperature of an ideal gas T^. The interaction 
strength is given by the scaling parameter rj = 0.31. Even though the scattering length is 
much smaller than the thermal wavelength, with a = 0.007At, the hard core repulsion sig- 
nificantly reduces the amount of bosonic attraction between the atoms. It becomes apparent 
from Fig. § that ^''•^-'(r) is dominated by gQ'\r){l — a/r)"^ corresponding to the first line of 
Eq. (|4.23|). The additional terms contribute a 1% correction. They mainly suppress the 



residual effect of the hard core repulsion for distances of the order of the thermal wavelength 
and larger. Fig. ^ also illustrates the diluting effect of the volume integration on the range of 
the bosonic attraction. For comparison, we plot in Fig. |]the correlation function g^"^^ for a 
Bose condensed gas at T = 0.5 Tq. Its maximum value is strongly reduced since condensate 
atoms do not exhibit any tendency to clump together due to wave function symmetrization. 
On the other hand, the existence of quasi-particle excitations leads to a relatively long tail 
for the correlation function. For this temperature, corrections to the term gQ'\r){l — a/r)"^ 



can be safely neglected, due to a larger thermal wavelength, i.e., a/Ar 
tail of the correlation function. 



0.005, and the long 



V. COHERENCE PROPERTIES 

We summarize in this section the coherence properties of trapped Bose gases based on our 
results for the first and second order correlation functions. First order coherence describes 
the ability of atoms with an initial separation r to form an interference pattern, when brought 
together. For an uncondensed gas, the visibility of such an interference pattern drops to zero 
on a length scale of r measured by the thermal wavelength, which is roughly of the order 
of 0.5/im in present experiments. In a Bose condensed gas, however, the visibility decreases 
on the same length scale to a plateau given by the total fraction of condensed atoms. This 
off-diagonal long range order allowed the experimental observation of an interference pattern 
I which extended over the whole cloud of condensate atoms. For spatial separations larger 



than the the size of the condensate, the visibility drops to zero as well, because long range 
order can not extend beyond the system size. 

Second order coherence implies instead the absence of any density correlations other 
than the trivial autocorrelation pointed out in Eq. (|2.8| ). In a dilute Bose gas, we can 
identify two sources of density correlations which limit the achievable degree of second order 
coherence. An inevitable anticorrelation between atoms originates from their strong mutual 
repulsion. The length scale of this repulsion is given by the scattering length a which is of 
the order of a few nanometers and therefore much smaller than any other relevant length 
scale. However, the relative wave function of two atoms increases as a/r, which is to say 
very slowly. Anticorrelation effects are therefore still visible on the length scale of the 
thermal wavelength, which is about 200 times larger (at T = 0.5 T°) than the scattering 
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length. The other hmitation to second order coherence is due to the tendency of bosonic 
atoms to cluster together. The range of this effective bosonic attraction is again limited 
by the thermal wavelength. However, its maximal amount is a function of temperature 
too, since the condensate atoms do not exhibit such a clustering effect. Because of the 
interplay between hard core repulsion and bosonic clustering, the maximum of the second 
order correlation function depends on both the temperature and the scattering length. 
It is nevertheless useful to characterize the local second order coherence by the maximum 

(2) 

of the fictitious correlation function (r, r) which does not account for hard core repulsion. 
We illustrate in Fig. ^ how this local measure of second order coherence, which is a function 
of the local condensate density, varies with the radial coordinate r in a spherically symmetric 
system (cf. ||18|). Analogously, the maximum 5'o^''(0) of the volume integrated correlation 
function offers a certain global measure of the second order coherence of the system. Fortu- 
nately, the quantity g^'' (0) is experimentally accessible due to its relation to the total atomic 
interaction energy. According to Eq. ( [4.18| ) the latter energy can be approximated by 

{H,)/2 = tr {h, /2 ^ tr {tG^ } /2 = ^J rfr G^^V, r) = ^ ^^^^(0). (5.1) 

The overlap integral between the true atomic interaction potential and the relative wave 
function between the atoms is therefore approximately proportional to the maximum value 
of the otherwise fictitious volume integrated correlation function Gq (r). 

In a recent experiment |jl5| the interaction energy and thereby gl)'\o) were determined 
by measuring the release energy of the atoms after switching off the magnetic trap. This 
experiment provided some evidence for a large degree of second order coherence of a trapped 
Bose-Einstein condensate at zero temperature. At finite temperature, however, a determi- 
nation of the interaction energy from the release energy would require that the initial kinetic 
energy of the uncondensed atoms be measured independently. In a similar experiment the 
third order analog of g^\o) has been determined by measuring three-body scattering rates 
16| . The implications of these experiments for second and third order coherence has been 



discussed before [15,16,18|- However, it was wrongly concluded that these experiments mea- 
sure the true second order correlation function G^'^^ or its third order analog. The significant 
impact of the hard core repulsion on second and third order coherence has apparently not 
yet been fully appreciated. 

Finally, we discuss the temperature dependence of g(^\o) in Fig. ^ The numerical plot 
shows a rapid but continuous decrease of (0) as the temperature falls below the crit- 
ical temperature. A comparison with a calculation in the Hartree-Fock limit shows that 
the existence of quasi-particle excitations has a small but noticeable effect on second order 
coherence. We contrast these results with the second order coherence of both a spatially ho- 
mogeneous and a harmonically trapped ideal Bose gas. In a spatially homogeneous medium 
one finds 









1 - 











(5.2) 



which is almost a linear function of temperature (cf. Fig. |^). In contrast the volume inte- 
grated second order correlation function of a harmonically trapped ideal Bose gas 
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/drG(i)(r,r)2' 
= l + 0(l/ViV), T<T, 



(5.3) 



experiences a sharp transition from the value two above the critical temperature to the 
value one below it. The reason for this discontinuous transition is that the peak density 
of a harmonically trapped noninteracting Bose gas becomes infinite as the thermodynamic 
limit is reached while the peak density of the thermal cloud retains a finite value. Of course, 
the presence of atomic interactions prevents any such unbounded growth of the condensate 
density for a real gas. Nevertheless, the spatial inhomogeneity accelerates the transition to 
a second order coherent ensemble as the temperature is decreased. 



VI. CONCLUSION 

We have extended the concept of optical coherence, defined in terms of correlation func- 
tions, to dilute trapped Bose gases. Except for a small region around the surface of a 
Bose-Einstein condensate, these correlation functions can be calculated by using a local 
density approximation. Even though the system behaves locally as if it were spatially ho- 
mogeneous, the spatial dependence of the effective chemical potential adds some interesting 
new physical features. The influence of both the Bose-Einstein distribution of uncondensed 
atoms and of quasi-particle like excitations is strongly reduced when a volume average over 
the mean coordinate of the correlation functions is performed. Nevertheless, consequences 
of quasi-particle excitations are still noticeable in the tails of the correlation functions. The 
density correlation is dominated by a tendency of bosonic atoms to cluster together and, on 
the other hand, the existence of a hard core repulsion between these atoms. The quantum 
statistical attraction and the effect of coulomb repulsion are found to be of comparable range 
for the conditions of most recent experiments. It is in this sense incorrect to neglect the 
hard core repulsion in the discussion of density correlations. 

For the discussion of second order coherence, however, it is convenient to introduce a 
fictitious correlation function which does omit the effect of hard core repulsion. This corre- 
lation function is experimentally accessible due to its close relation to the atomic interaction 
energy. By introducing this fictitious correlation function, we are able to explain the re- 
sults of recent experiments on second and third order coherence. Finally, we have observed 
that a harmonically trapped gas achieves second order coherence much more rapidly as the 
temperature is lowered below the critical temperature than does its spatially homogeneous 
counterpart. 
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FIG. 1. This plot shows the volume inte- 
grated first order (bottom) and second order 
(top) correlation functions for a harmonically 
trapped ideal gas at the critical temperature. 
The solid lines correspond to the exact expres- 
sions of Eqs. ( p2D and (|^, while the dashed 
lines approximate the occupation distribution of 
the trap levels by a Maxwell-Boltzmann distri- 
bution as in Eq. ( 3.23| ). In order to display the 
diluting effect of the volume integration, implicit 
in these results, we have added as dot-dashed 
lines the unintegrated correlation functions for 
mean positions at the center of the trap, where 
the local chemical potential is the smallest. 



FIG. 2. The solid line shows the volume inte- 
grated first order correlation function for a har- 
monically trapped interacting gas at T = 0.5 
where is the critical temperature of the ideal 
gas. The scaling parameter is 77 = 0.31. The 
correlation function does not decay to zero due 
to the existence of a Bose-Einstein condensate. 
As in Fig. ^, the dashed line corresponds to a 
Maxwell-Boltzmann distribution of the uncon- 
densed atoms while the dot-dashed line gives the 
unintegrated correlation function at the center of 
the trap. The existence of quasi-particle excita- 
tions leads to a significantly slower decay of the 
tail of the correlation function, compared to the 
case of an ideal gas. 
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FIG. 3. The volume integrated second order 
correlation function of an interacting gas is given 
by the solid line. The temperature corresponds 
to the critical temperature of an ideal gas. 
With a scaling parameter of r/ = 0.31, the scat- 
tering length a is much smaller than the ther- 
mal wavelength, a = 0.007 At. Nevertheless, the 
hard core repulsion leads to noticeable reduction 
of the correlation function over a broad range of 
relative distances. The dashed line shows as a 
comparison the volume integrated fictitious cor- 

(2) 

relation function (r) which omits the effect 
of the hard core repulsion. 
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FIG. 5. The volume integrated second order 
correlation function is here plotted for a temper- 
ature T = 0.5 T^. The dashed line shows gl^\r) 
which omits the hard core repulsion. Due to 
the presence of a considerable condensate frac- 
tion, a large degree of second order coherence is 
attained. The existence of quasi-particle excita- 
tions leads to a remarkably slow decrease of the 
correlation function. 
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FIG. 4. This plot compares the long range 
behavior of the volume integrated second order 
correlation function of Fig. ^ (solid line) with its 
main constituent g^\r){l — a/r)"^ (dashed line). 
The latter approximates the relative wave func- 
tion by its zero temperature limit. The differ- 
ence is small, and becomes even smaller at lower 
temperatures. The dot-dashed line, by contrast, 
displays the unintegrated correlation function at 
the center of the trap. 
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FIG. 6. The normalized second order corre- 
lation function (r, r) of Eq. (3.34) is a mea- 
sure of local second order coherence. It also 
describes the ratio of the total interaction en- 
ergy, including quantum statistical exchange en- 
ergy, to the classical Hartree mean-field energy. 
We have plotted (r,r) (solid line) for a har- 
monically trapped gas with interaction strength 
r\ = 0.31 and temperature T = 0.8 T^. For com- 
parison, the dot-dashed curve shows the total 
density distribution in arbitrary units. 
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FIG. 7. The peak value g^ (0) of the vol- 
ume integrated second order correlation func- 
tion, excluding the hard core repulsion between 
the atoms, is a good measure of the second or- 
der coherence of the system. The solid line shows 
the influence of the temperature on the second 
order coherence for a harmonically trapped in- 
teracting gas in the Popov approximation with 
scaling parameter r] = 0.31. A corresponding 
calculation using the Hartree-Fock approxima- 
tion (dashed line) shows that quasi-particle ex- 
citations have a small but noticeable effect. For 
comparison, the dotted line describes a spatially 
homogeneous ideal gas while the dot-dashed line 
corresponds to a harmonically trapped ideal gas. 



25 



